Scaling laws for weakly interacting cosmic (super) string and p-brane networks 
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In this paper we find new scaling laws for the evolution of p-brane networks in A'^ + 1-dimensional 
Friedmann- Robertson- Walker universes in the weakly-interacting limit, giving particular emphasis 
to the case of cosmic superstrings (p = 1) living in a universe with three spatial dimensions (A^ = 3). 
In particular, we show that, during the radiation era, the root-mean-square velocity is w = 1/^/2 
and the characteristic length of non-interacting cosmic string networks scales as L oc a^^^ (a is the 
scale factor), thus leading to string domination even when gravitational backreaction is taken into 
account. We demonstrate, however, that a small non-vanishing constant loop chopping efficiency 
parameter c leads to a linear scaling solution with constant LH <^ 1 [H is the Hubble parameter) 
and u ~ 1/ y/2 in the radiation era, which may allow for a cosmologically relevant cosmic string role 
even in the case of light strings. We also determine the impact that the radiation-matter transition 
has on the dynamics of weakly interacting cosmic superstring networks. 



I. INTRODUCTION 

Cosmic strings are usually regarded as benign objects 
with potentially interesting cosmological consequences. 
Ordinary cosmic strings have an efficient energy-loss 
mechanism: once two strings collide, they exchange part- 
ners and intercommute. This results in the creation of 
string loops, which eventually decay radiatively. The cor- 
responding energy loss aids the attainment of a linear 
scaling regime, during which the average energy density 
of the string network remains a constant fraction of the 
background energy density . The existence of this 
regime has not only been confirmed numerically [3|-l5|, 
but was also established using semi-analytical models @- 

Recent developments in string theory indicate that fun- 
damental strings (F-strings) or 1-dimensional Dirichlet 
branes (D-strings) may grow to macroscopic scales and 
play the role of cosmic strings p^l - [l3l |. These cosmic su- 
perstrings may have an intercommuting probability sig- 
nificantly smaller than unity [T^ . and thus cosmic su- 
perstring networks may be weakly interacting. It is then 
crucial to understand the effect that a significantly less 
efficient energy-loss mechanism might have on the late- 
time evolution of the network, and to determine the con- 
ditions under which it might prevent the string network 
from dominating the energy density of the universe. 

In this article, we revisit the question of whether or not 
weakly interacting string networks are able to attain lin- 
ear scaling regimes. We start by reviewing, in Sec. |lll the 
dynamics of p-branes and the Velocity-Dependent One- 
Scale (VOS) model for p-brane networks. In Sec. IIII[ we 
obtain the frictionless scaling laws for weakly interact- 
ing p-brane networks. In Sec. llVi we consider the par- 
ticular case of cosmic superstrings in 3 -t- 1-dimensional 



Friedmann-Robertson- Walker (FRW) backgrounds, and 
determine the impact that gravitational backreaction has 
on the dynamics of the network. We then conclude in Sec 
lYl 



II. p-BRANE NETWORK DYNAMICS 

In the zero-thickness limit, the world- history in space- 
time of a featureless p-brane may be represented by 



(1) 



where u" with = 0, 1, ...,p are the coordinates parame- 
terizing the (p -|- l)-dimensional worldsheet swept by the 
p-brane, vP is a timelike parameter and are spacelike 
parameters. The action of a thin and featureless p-brane 
is described by the Nambu-Goto action 



S = -Or. 
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where Up is the (constant) p-brane mass per unit p- 

dimensional area, g = det(g^c;), g^jj — gajix^jj^x^- is the 

worldsheet metric, x'^ — dx'^ jdv?- and g^j^i, is the metric 
tensor. 

In a {N 4- l)-dimensional flat Friedmann- Robertson- 
Walker Universe, the line element is given by 



a^irj) [d^rj — dx ■ dx) , 



(3) 



where a represents the cosmological scale factor, t and 
rj = J dt/a are respectively the physical and conformal 
time and x is a A^-vector whose components are comov- 
ing cartesian coordinates. By varying the action in Eq. 
^ with respect to and imposing temporal-transverse 
gauge conditions. 
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one obtains the equation of motion [15 1 



VOS model 



P r 
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Here 



1 



^2 



(6) 
(7) 

(8) 



a dot represents a derivative with respect to conformal 
time, = dx/du^, are orthogonal spacehke world- 
sheet coordinates (with i = I,-- - and H — a/ a. 
Given the gauge conditions described in Eq. (O, x rep- 
resents the physical velocity of the p-brane and it is per- 
pendicular to the brane itself. 

The total energy E and root-mean-square (RMS) ve- 
locity t; of a p-brane network are defined respectively by 



E 



JedPu ' 



so that 



E = {p + i)nE ( 



p 



\p+l 



(9) 
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For very small p-branes it is a good approximation to 
consider that the expansion of the universe has, on aver- 
age, no impact on the total energy, so that the averages 
over a sufficiently long time of the total energy {E)t and 
RMS velocity {v'^)t = p/{p+ 1) are constant. 

In the case of a statistically homogeneous p-brane net- 
work, the average p-brane energy density may be defined 
as pp = E/V, with V cx a^. It is then straightforward 
to show that 



Pp + npp [D+{p+ 1) v^]=0, 



(12) 



where D — N — p. By defining the characteristic length, 
L, of the network as 



Pp - 



this equation may be rewritten as 
dL 



where H — T-L/a is the Hubble parameter. 



(13) 



(14) 



A unified VOS model for the dynamics of p-brane net- 
works in (iV+ l)-dimensional FRW universes was derived 
in [i^ According to this model, the cosmological 

evolution of a statistically homogeneous p-brane network 
may be described by the following equations 



au , 



dv 
~dt 
dL 

dt 



V 



= HL 



L 
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(15) 
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where £^ ^ {p + \)H is the damping lengthscale, fc is a 
dimensionless curvature parameter and the assumption 
that (v'^) — v"^ was made (see §j). Eq. ([T5| is iden- 
tical to Eq. ([T^ apart from the energy loss term asso- 
ciated with p-brane reconnection (last term) while Eq. 
(ITSI) has been rigorously derived using Eqs. ([5]), ^ and 
(flUl) in [T5l - [l7| (the curvature parameter k is also rigor- 
ously defined therein). A frictional force — caused by 
the interaction of the branes with ultrarelativistic parti- 
cles or other frictional sources — may also be included in 
Eq. ((T5|) , by introducing a friction lengthscale so that 



{p + l)H 



For simplicity, we will not con- 



sider it in the present paper and consequently we shall 
take £f — oo (see [13] for a detailed discussion of friction 
dominated regimes). 



III. SCALING REGIMES 

In a FRW universe with a decelerating power-law ex- 
pansion — aoct^ with < /3 < 1 ~ Eqs. ([151) and (fT6)) 
admit linear attractor solutions of the form 



and 



V = constant , 



with 



k{k + c) 



/3{l-(3)D{p + l) 

I (1 - l3)kD 
I3{k + d){p+l) ■ 



(17) 

(18) 
(19) 



The existence of an energy-loss mechanism helps the at- 
tainment of a linear scaling solution but that is not a 
necessary condition. It is a common misconception that, 
if c = 0, these solutions are possible for all /3 > _D/(A^-f-l) 
(or, equivalently < 1 in Eq. (|19p ). However, in ex- 
panding backgrounds, it is reasonable to expect the RMS 
velocity to be smaller than it would be in Minkowski 
space. One should then expect k{v) to vanish for — 
p/{p+l) and the scahng velocity not to exceed this value. 
This, in turn, implies that the linear scaling solutions are, 
in the absence of energy-loss mechanisms due to p-brane 
reconnection, attainable only for (3 > D/N. In a 3 -I- 1 



3 



dimensional FRW universe, the threshold is then /3 = 2/3 
(rather than /3 — 1/2) for cosmic strings, and /3 — 1/3 
(instead of /3 = 1/4) for domain walls. 



may be found up to first order in c and /. Note that, in 
the 13 D/N limit, one needs to use Eqs. (IT71)-((TH]) in 
combination with Eqs. (|23p and (|24p in order to find the 
scaling solutions. 



A. Non-interacting p-branes 

For c — and < (3 < D/N, p-brane networks are 
expected to evolve towards a solution with ~ p/{p+l). 
The equation-of-state parameter of a brane gas is given 
by 



Pp 



1 

N 



[{p +l)v^- p] 



(20) 



where Vp is the average brane pressure. One may then 
conclude that, for < /3 < D/N (and = p/{p+l)), the 
network can be described, on scales much larger than L, 
as a quasi-homogeneous matter background with Wp = 0. 
If c = 0, Eq. (|16p may be written as 



(21) 



We then have that, \iQ < (3 < D/N (and = p/(p + 
1)), the characteristic lengthscale of the network scales as 
L (X a^^^ . Hence, in this case, the characteristic length 
evolves asymptotically as 



with 



(22) 



B. Weakly interacting case 



The picture changes significantly if we consider a con- 
stant non- vanishing energy loss parameter < c ^ 1. 
This case is relevant not only for the study of cosmic su- 
perstring networks — which might have a small intercom- 
muting probability — but also for the study of p-brane 
networks in iV-dimensional backgrounds. In particular, 
if p < {N — l)/2, p-branes are likely to miss each other 

m- 

As previously discussed, the curvature parameter 
should vanish and change sign at v'^ = p/{p + 1). Let 
us define the parameter k(v) as 



k{v) = Af{v) , 



where A is a constant and 



fiv) 



p+1 



(23) 



(24) 



The curvature parameter defined Eqs. ([25)) and (IM|) may 
be considered a good approximation if / ^ 1. In this 
case, for c ^ 1 and /? < D/N, a linear scaling solution 
(see Eq. ^) with 



(p+ 1)1/2(1) _ /37V) 



fl„l/2(„+ 1)1/2 
/ = Y ^ ^^^^ 



IV. COSMIC STRINGS 

Let us consider the more interesting case of cosmic 
superstring networks in a 3-1-1 dimensional FRW back- 
ground (with N = 3, p = 1, and /i = cti). The generaliza- 
tion to other values of p and N is trivial. It follows from 
Eq. (j22p that during the radiation era, with /3 = 1/2, 
one has L oc c?^"^ oc i'^/'*, in the absence of any energy- 
loss mechanisms. Therefore, the string energy density 
Pfi — Pi — p/L^ is proportional to and tends to 
dominate the energy density of the universe during the 
radiation era. 



A. Gravitational Backreaction 

The effect of gravitational backreaction on the net- 
works' dynamics may be taken into account by adding 
the following term to the right-hand side of Eq. (TTF 



dL 



(26) 



with r ^ 65 [19l - l21j . In this case, although a scaling 
solution given by Eq. (jl7p with 



would be mathematically correct, it would imply 



Pb / rad 



(27) 



(28) 



where pt is the background density during the radiation 
dominated era {pb = l/(327rGt2)). However, this possi- 
bility is completely excluded by observations: the cosmic 
string tension is constrained to be Gp < 10~^ p3 - [2^ . 
Thus, for any reasonable value of Gp, the right hand 
side of Eq. (j28p is much larger than unity. The inclusion 
of the effects of gravitational backreaction does not pre- 
vent cosmic strings from dominating the energy density 
of the universe. 

On the other hand, one might consider a constant non- 
zero energy-loss parameter such that Gp ^ c <C 1 and 
take the momentum parameter proposed in 9] 



TT 1 



8v^ 



(29) 



Then, during the radiation era, the VOS equations admit 
a linear scaling solution as in Eq. (jl7p . with 



C = V2c 



1 

71 



(30) 
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and 



c 



(31) 



up to first order in c (here B = 12/tt). Note that during 
this regime v is not exactly 1/V2. The deviation 6 was 
be estimated by proceeding in a similar fashion as in Sec. 
IIII Bl and by noting that, up to first order in (1 — v^w), 
Eq. ([^ may be written as 



(32) 



k{v) = B[l/V2-v 



During this regime, one has 



167r 



G/i 



(33) 



rad 



which means that cosmic strings could make a fairly large 
contribution to a non-clustering matter background dur- 
ing the radiation era, even for low values of G/z, as long 
as c is small enough. For instance, if c = 10~^ and G/i = 
10"-^^ then (/Op/pfc)j.ad = 5 x 10^''. On the other hand, 
if g = IQ-^ and G/i'^ 10"^ then (p^/pfc),^^ = ^ x 10"^ 
Note that c is expected to be proportional to the string 
intercommuting probability which may be significantly 
smaller than unity. According to values of c within 
the range [10~*, 10~^] appear to be acceptable (c ~ 0.23 
if the string intercommuting probability is equal to unity 
i). 



B. Radiation-Matter Transition 

Let us consider a fiat FRW universe containing only 
matter and radiation, thus neglecting the dynamical ef- 
fects associated with the recent acceleration of the expan- 
sion of the universe. Using the VOS equations (Eqs. 
and (jl6p ). it is possible to find numerical solutions for 
the evolution of non-interacting and weakly interacting 
cosmic string networks throughout the radiation-matter 
transition. The results are illustrated in Figs. [T] and [2j 
which show the evolution oi ^ = L/t and v for a network 
of cosmic strings in a fiat 3 + 1-dimensional FRW uni- 
verse with Gfj, = 10~^^ and c = (solid (purple) line) or 
c = 10~^ (dash-dotted (blue) line) as a function of the 
scale factor, a (the value of the scale factor at radiation- 
matter equality is normalized to unity). 

Figs. [T] and m show that, as previously discussed, cos- 
mic string networks are unable to attain a linear scal- 
ing regime during the radiation era if c = 0. Instead, 
the RMS velocity and the characteristic length of the 
network evolve rather quickly towards a scaling regime 
with ■u^ = 1/2 and ^ = L/t (x i"^/^ cx a'^/^. During 
the radiation-matter transition the rate of change of ^ 
decreases and, eventually, tends to a constant. The 
network attains a linear scaling regime, with constant 
^, only when the universe becomes matter-dominated. 
Note, however, that the scaling value of ^ is not defined 




a = . Gii = iQ-'-i- 
c = 10 ^. G/i = 10" -• 



FIG. 1: Evolution oi ^ = L/t for a network of cosmic strings 
in a flat 3 + 1-dimensional FRW universe with Gjj, — 10~^^ 
and c = (solid (purple) line) or c = 10~^ (dash-dotted (blue) 
line) as a function of the scale factor , a. The (black) dotted 
horizontal line indicates the expected value of ^ during the 
radiation era (see Eq. (|30|) '). The value of the scale factor at 
radiation-matter equality is normalized to unity. 
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FIG. 2: Evolution of the RMS velocity, v, for a network of 
cosmic strings in a flat 3 + 1-dimensional FRW universe with 
Gfj, = 10~^^ and 5 = (solid (purple) line) or c = 10~^ (dash- 
dotted (blue) line) as a function of the scale factor, a. The 
(black) dotted horizontal line corresponds to i; = l/\/2. 



in this case: it is simply the value of ^ at the onset of 
matter domination and consequently it depends on the 
initial conditions. 

On the other hand, if energy loss due to loop forma- 
tion is taken into consideration (with G/i ^ c <C 1), then 
cosmic string networks may undergo two separate linear 
scaling regimes. If the string-forming phase transition 
occurs deep enough in the radiation era, the network is 
able to attain the regime in Eq. ((30|) . As Fig. [2] clearly 
illustrates the value of v is fairly close to w ~ l/v^ in 
the radiation era, but there is a slight deviation. In this 
particular case, the value of S is in good the agreement 



5 



with the theoretical prediction in Eq. (PT|) : the rela- 
tive difference is smaller than 1%. As the universe en- 
ters the matter dominated era, the network undergoes a 
phase during which ^ grows (quasi-logarithmically) and 
V decreases as time progresses, until it reaches the linear 
scaling solution defined in Eqs. (fT7|) - ([T^ with k{v) given 
in Eq. ([29]) . Note that, the smaller the value of c, the 
longer it takes the network to attain the scale invariant 
regime in the matter-dominated era. For weakly inter- 
acting networks the transition phase may indeed be long 
lasting. 

V. CONCLUSIONS 

In this paper we determined the macroscopic evolution 
of p-brane networks in iV -I- 1-dimensional FRW universes 
with a decelerating power-law expansion a oc {0 < /3 < 
1). We found that, for /3 < D/N, a linear scaling solution 
with L oc t is possible only if the energy loss term due 
to p-brane reconnection and decay (parametrized by c) 
is non-zero. We have shown that, in this case, if c = 
then L cx a^^-^ oc t^^^^ , eventually leading to a p-brane 
dominated universe. We have further shown that if /3 < 
D/N and the p-brane network is weakly interacting (with 
a small but non-zero c) the network then evolves towards 
a linear scaling solution with ?; ~ p/(p-|- 1) and LH <C 1. 
The corresponding equation-of-state parameter is Wp ^ 
thus generating a quasi-homogeneous matter background 



on scales much larger than L. 

Our results have profound implications for the dynam- 
ics of cosmic superstring networks in expanding 3 + 1- 
dimensional FRW universes. If c = then the root-mcan- 
square velocity tends to v = 1/V2 and the characteristic 
scale of non-interacting cosmic string networks scales as 
L (X a'^/^ during the radiation era. We have shown that 
this leads to string domination even when gravitational 
backreaction is taken into account. On the other hand, 
a small non-zero constant loop chopping efficiency pa- 
rameter c leads to a linear scaling solution with constant 
LH <^ 1 {H is the Hubble parameter) and v ~ l/\/2 in 
the radiation era, thus generating a quasi-homogeneous 
matter background on scales much larger than the char- 
acteristic length L of the network. We have shown that, 
in the scaling regime, the ratio between the average string 
energy density and the total background density is pro- 
portional to G/i/c" which may allow for a cosmologically 
relevant cosmic string role even for light strings. We have 
also found that the transition from the radiation to the 
matter dominated era has a major impact on the dynam- 
ics of weakly interacting cosmic superstring networks. 
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